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We consider a model of the dynamics of a magnetic 
grain, incorporating interactions with nuclear and paramag- 
netic spins, conduction electrons, and phonons. Decoherence 
comes both from the spins and the electrons, although elec- 
tron effects can be neglected if the substrate is insulating, and 
the grain has spin S < 10"". A successful experiment will re- 
quire both isotopic and chemical purification, if coherence is 
to be seen; we suggest a "spin echo" arrangement to do this. 

PACS numbers: 75.10.Jm, 75.60.Jp 

Much activity in the new field of nanomagnetism |^] 
focuses on the quantum dynamics of small magnetic par- 
ticles or wires [|[||. Of potential interest to the mag- 
netic recording, information storage, and computer in- 
dustries, this work also seeks to understand quantum 
magnetism on the mesoscopic scale - practical applica- 
tions are still unclear. Recently a debate has arisen over 
experiments claiming to see quantum coherent mo- 
tion of the magnetisation in ferritin nanomolecules. Criti- 
cisms have centered around, e.g., (i) the dipolar coupling 
between molecules (ii) the experimental power ab- 
sorption [Hj^ , (iii) decoherence or coherence blocking by 
environmental spins [p|~p^ , and (iv) the apparent incon- 
sistency with high-T magnetic blocking |jll|. Neverthe- 
less recent experiments ||^ see further evidence of coher- 
ence. 

Here we present a theory of coherence in mesoscopic 
magnetic molecules or grains, including environmental 
spins, phonons, and possible electrons, and also the cou- 
pling to the measuring system. This has implications for 
the ferritin experiments; more importantly, it shows what 
has to be done in future experiments to see coherence. 

Many of our results can be understood qualitatively 
by simple arguments. At low T, one models the grain 
as a "giant spin" S, with effective Hamiltonian trun- 
cated to the 2-level 11° fj — 2Aofj: costtS, involving co- 
herent motion between, e.g., the 2 quasiclassical states 
S — ±zS at a frequency 4Ao when symmetry permits 
To maintain coherence one requires (i) near de- 
generacy, within energy ^ Ao, of the 2 states, and (ii) no 
phase decoherence. 

Consider now a Hamiltonian 
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with hyperfine couplings {uJk} between S and nuclear 
spins {/fc}; with no nuclear spins Aq = r^oC^'^", where 
VLo - 2{K\\K^Y/^, and Ao - S{K\\/ Ka_Y/^ . Typically 
^o ~ 10^ - 10" Hz, and ~ 10^ - 10^ Hz; for meso- 
scopic spins [S > lO'^), one expects Aq < 10^ Hz, for 
tunneling between ztzS", so Aq <C LOk- 

However, even if the applied field Ho = 0, an inter- 
nal bias field e — Y^^^i^klk acts on S. Suppose we 
have a dominant hyperfine coupling Wq, with a spread 
Sujk of levels around lOo, caused by Nakamara-Suhl in- 
teractions {6u)k ~ 10"^ — 10^ Hz), transfer hyperfine cou- 
plings {SiJk ~ 1 — 20 MHz) or other magnetic hyper- 
fine couplings. Define a "density of states" W{e) for the 
bias; if /i = N^^^SuJk/(^o ^ 1, this consists of "polarisa- 
tion groups" of width ~ /iWo around e = lUoAN, where 
AiV — — is the total nuclear polarization [|l^, in- 
side a Gaussian envelope of width N^^^lUo. Usually how- 
ever 3> 1, and different polarization groups overlap. 
Notice only the fraction A of grains in an ensemble hav- 
ing e < Ao can flip - the rest are " degeneracy blocked" . 
Thus if ^ <C 1, a fraction Ao/{N^/^SLJk) of that por- 
tion / = ^/2/^ of grains having AiV = can flip, so 
A - Ao/{NSujk); ff A* > 1, then A ~ Ao/iN^/^tOo). 

Even these few grains may not necessarily flip co- 
herently. Flipping S causes transitions in the nuclear 
system, which destroy phase coherence, even if the 
change in polarisation AM — 0; and other fields can 
cause a mismatch between initial and final state nuclear 
wave-functions, analogous to the Anderson orthogonal- 
ity catastrophe |14|. Electrons or phonons can remove 
the degeneracy blocking, by absorbing the bias energy, 
but further destroy coherence in doing so. Finally, nu- 
clear spin diffusion, at a rate Tj"^, destroys coherence by 
causing e to fluctuate in time. 

To see coherence is clearly going to be hard; we now 
proceed to a quantitative theory, to see what can be done. 

(i) The Spin Environment: At energies <^ flo, spin ef- 
fects are treated by truncating to an effective Hamilto- 
nian Heff = HsN + Hnj\j, where [|8|-[l0|] 
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HsN describes all couplings between S (now r) and the 
spin-1/2 environmental variables {^k}- The dimension- 
less parameters ak, ^k, and <i>, the energies Ao, wf, and 
uj^, and the unit vectors fik, Vk, Ik, iTik, have been de- 
rived for various "bare Hamiltonians" . If, e.g., the {Ik} 
in d^) are spin-1/2 nuclei, one finds Aq = Aq, $ =TrS, 
= ak/V^ = 7ra;fc/20o for small cj^/fio (for general 
Wfc/rio, see Refs. |||,|l^), w| = lu^, and uj^ = 0, whilst 
fik — 2~^/^(x, y), Vk = —X, Ik — z, and mfc is arbitrary. 
The "longitudinal coupling" wjl gives the change in en- 
ergy of (Tfc when S flips, and to'^rhk represents any field on 
(Tfc, perpendicular to Ik, which is unchanged by tunnel- 
ing; this arises if ak comes from truncating a higher-spin 
system (e.g., if Ik is spin 1), or if the 2 quasi-classical 
orientations of S are not exactly antiparallel. In reality 
Q!fc <ti 1, and is just the amplitude for ak to coflip with 
5; when S flips, roughly A = 1/2 J^k '^1 nuclei also flip. 

To study coherence we calculate P{t) = {Tz{t)fz{0)) 
[ p5[ , which shows coherent oscillations without the envi- 
ronment (i.e., P{t) P(°)(t) = cos[4Aoicos7r5]). We 
have solved for P{t), including all terms in (|^) and <^). 
We write it as 
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with W{e) as before, and Z{T) the nuclear partition func- 
tion. Eq.(^ averages over bias e and sums over pro- 
cesses in which the polarisation changes by 2M . The 
case uj^ = ^fe = is sufficient M] for what follows; then 
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where if ^ = (no spin diffusion) , we have 
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describing a biased 2-level system, where E^j = e 



and AM{'P,y) = 2Ao | cos (/?Jm (2^ (A - X')y) |, with 
A' = l/SEfcafel'^fc)^ (so that A > A'). The averages 
over topological phase ip, winding number v (cf. [g[), 
and "orthogonality blocking" y, then describe all the dy- 
namical effects of the {<Tk} on S. The weighting fac- 
tor Fy{i/) = exp{-4AV^}; often, as for (§), nl = 0, so 
A' = and the phase average collapses to a delta- function 

In rare cases A ^ 1 (i.e., no nuclei flip), and (still 
assuming Tj"^ = 0) the nuclear bath acts as a static field. 



An example may be ferritin, where only 2.16% of the Fe 
nuclei have spin-1/2 (so N ^ 100), and uJo ~ 50 MHz 
(but fi > 1 if we take account of hyperfine interaction 
with 5000 nuclei). Assuming Q that fio 4 x 
10^°Hz, one gets A - 4 x 10"^ Then ff ksT > lUo, 
Eq.(||) collapses to P{t) = I - ^AJ^k J2k+i['iKt], with 
A = 7rAo/(woV27riV) (for ferritin, a A^ ~ 1 MHz is 
claimed, so ^ ~ 2 x 10"'^), equivalent to an absorption 
spectrum (Fig.l): 
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i.e., the spectrum of a set of biased 2-level systems with 
a fraction A near resonance [T^ ] . There is no dissipation, 
but the spectrum is not a sharp line . 



FIG. 1. The frequency response x"{^) for the giant spin 
assuming Tj"^ = 0, for the cases (a) when A = 10~^, and only 
degeneracy blocking occurs, and (b) when A — A' = 10, and 
topological decoherence also plays a role. 

Unfortunately in reality T^^ ~ 10^ - 10*^ Hz (the 
strength of the Nakamura-Suhl interaction); this causes e 
to fluctuate in time, with ((e(t) - e(0))2)V2 _ utOoT^h, 
over the energy range ficuo of the polarisation group in 
which the nuclei lie. Thus in a time A^^, e{t) changes 
by Ae ~ iiujo{T2Ao)~^^^ ■ Coherence will be destroyed 
unless Ae ^ Aq, i.e., coherence requires 



(8) 



where the latter expression assumes ~ Sujo, i.e., that 
SuJo is mostly dipolar spreading [jisl . Thus it is crucial 
that the spin diffusion rate Tj"^ be kept well below Aq. 
In ferritin, most nuclear spin diffusion will occur amongst 
the N ~ 5000 protons; from (H) we see that coherence 
then requires <C 10^ Hz, otherwise the peak in (|^) 
will be entirely washed out (we are unaware of any T2 
measurements in ferritin). 

In fact for most grains A ~ 0(1) or greater, and the 
coherence is destroyed simply by the coflipping of nuclear 
spins (i.e., by topological decoherence |^). Fig.l shows 
the case A = 10, appropriate to Tb grains containing 
- 10^ Tb ions. 

So much for nuclear spins inside the grain - what of 
paramagnetic impurities and nuclei in the surrounding 
substrate, which feel the dipolar field of the grain? These 
dipolar fields are small compared to fio, but for a meso- 
scopic grain, their number N will be very large, and if A 
is ~ 0(1) or greater, then they too will destroy coherence 

(ii) Phonons fc Electrons : In the presence of environ- 
mental spins, electrons and phonons only further hinder 
coherence. Suppose however we can eliminate nearly all 
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such spins. Consider now the electrons - these are partic- 
ularly dangerous for coherence. If grain and substrate are 
conducting a reasonably good description of the coupling 
to S is 
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where the form factor — J^{d^f/Vo e*^'"p(r) inte- 
grates the number density p{r) of grain spins over the 
grain volume VJ, ~ -Ri?, and ct creates electron mo- 

kjCt 

mentum and spin eigenstates; by standard techniques 
we can derive a Caldeira-Leggett spectral func- 
tion J{uj) = TTabio, of Ohmic form, for this coupling, 
where at = ^g^S^d F^^p \^)f.s. averages | p over 

the Fermi surface (| k |, | k' |~ kp), and the dimension- 
less Kondo coupling g = JN{0) ^ 0.1 for metals; then 
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SO ab ^ 1 for mesoscopic grains [£l|, completely de- 
stroying coherence at any temperature [ jist . Even if the 
grain is insulating, the surface spins still couple to sub- 
strate electrons, and ( |To| ) is replaced by a surface coupling 
as g^S^^^, with probably smaller exchange coupling 
gs- Coherence then requires the damping T ^ ~ 2'KaskBT 
to be considerably less than the renormalised splitting 
Ar = Ao(Ao/rio)"''^^"'^~"''\ which is hard to imagine for 
any reasonable ksT unless S* ^ 100. 

The case of conducting grains on an insulating sub- 
strate is more interesting; now the electronic states are 
discrete, with separation Ae ^ D/S > 10"^ /S in tem- 
perature units, where D is the bandwidth. The de- 
tailed theory, involving both mesoscopic level fluctua- 
tions and surface scattering is complicated, but the re- 
sults are intuitively obvious. If S flips, the electronic 
spins follow adiabatically since J 3> 5^o, but surface spins 
scatter electrons from one internal orbital state to an- 
other. This gives an Ohmic coupling with ~ Q^s/[1 + 
exp{Ae/kBT)], and a renormalized splitting A^ ~ Aq (if 
Ae > no), or Ar - Ao(Ae/rJo)"% if Ae < ^o- The 
damping rate is now ~ 2-K'askBT , and rewriting Fj, 
Ar, and Ae in terms of 5*, one finds almost perfect co- 
herence provided that 
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For experiments at mK temperatures, this means that 
electronic dissipation effects can be ignored in such grains 
provided S" < 10^ 

For phonons the damping rate F^ is p^j20[ | 



T^{kBT)^(^^ J(A<,) coth(A„/2fcBT) , (12) 

where J{uj) - Ciu/Qof] the coupling C ~ [KSY/Qd, 
where K is a, typical anisotropy energy, and Od is the 
Debye temperature. For the relevant magnetoelastic in- 
teractions the effect of phonons on tunneling from a 
metastable well was investigated in ppl , but here there is 



an extra factor (Ao/Slo)^ which arises |20 because there 
is no "diagonal coupling", i.e., no coupling to our from 
phonons (in zero external field) . Even without this factor 
F^ is very small (inverse minutes or hours). There will 
also be Ohmic 2-phonon couplings, but these give pC| ] 
an even smaller damping ^ [kBT /QoY ■ Thus, as usual, 
phonon effects on coherence are negligible |^3j . 

(iii) Real Experiments: We now see that to see coher- 
ence one requires, at the least (i) isotopic purification 
of the magnetic grains, and (ii) an insulating substrate, 
purified of paramagnetic impurities. The grains may be 
conducting if they are not too large - in any case they 
should be purified of paramagnetic impurities. 

These sample preparation requirements do not seem 
to be beyond current techniques. This only leaves de- 
coherence coming from the measuring apparatus, which 
also interacts with the grains. One solution has been 
discussed by Tesche p^ , for coherence experiments in 
squid's. Here we propose another "spin-echo" set-up, 
which might be easier to implement for magnetic grains. 
If there is a sharp line in the absorption spectrum due 
to coherence, then a long-lived echo signal should exist. 
Moreover, even if the grains have a wide distribution of 
bare tunneling rates Aq, so x"(w) is basically structure- 
less, spin echo and "burned hole" type experiments may 
still unveil the existence of coherence ||2^. For example, 
by applying the standard 7r/2 magnetic pulse (to gener- 
ate nondiagonal elements of the giant spin density ma- 
trix) followed by one or more tt and 7r/2 pulses to reverse 
the phase decoherence due to inhomogeneous broaden- 
ing and/or introduce a time delay in the transverse re- 
laxation, one can extract the magnetic moment of the 
grains and check whether it is consistent with suscep- 
tibility experiments on the same grains. One may also 
study damping mechanisms by observing the shape and 
amplitude of the "burned hole" or multiple echo signal. 
These depend, of course, on the distributions of Aq and 
e and on the decoherence mechanism - the details will be 
given in a longer paper [Q. 

The crucial points here are (1) coherence can be found 
even for a spread of {Aq} over grains, (2) spurious peaks 
in x" can be eliminated, and (3) we have a non-invasive 
measurement, where the operator exp{— j J H(T)dT}, 
acting between the pulses, does all the work. This last 
point is of course very important within the context of 
quantum measurement theory. 
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